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Abstract 

We analyze the nested derivatives of a function 2) n [/] (x) asymp- 
totically, as n — > oo, using a discrete version of the ray method. We 
give some examples showing the accuracy of our formulas. 
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1 Introduction 

The oldest and most widely used method for computing the Taylor series of 
inverse functions is the Lagrange Inversion theorem [TTJ , which can be stated 
as follows [TT] : 

Theorem 1 Suppose that f(z) is analytic at a, and f'(a) ^ 0. Then, 

oo 

f- 1 ( z ) = a + J^c„ (z - b) n , 

71=1 
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on a neighborhood of b = f(a), where 



1 d n 1 w — a I™ 
n n\ dw"- -1 f(w) — b 

Modifications and extensions of this formula were studied by Apostol pQ , 
Gessel [SJ, Krattenthaler [TO], Roman and Rota [15] and Sokal [T6j among 
others. 

Formulas relating the coefficients of the Taylor series of a function and 
its inverse were obtained by Jacobsthal [9], Rauch [H] and Ostrowski [T2] . 
The two- variable case was considered in |13j . 

In [2], we derived an algorithmic approach that simplifies the computation 
of the derivatives of inverse functions. We defined 2) n [/] (x), the n th nested 
derivative of the function f(x), by D°[f] (x) = 1 and 

2> n+1 [/] (x) = ^ [f(x)V n [f] (x)} , n = 0,l,.... (1) 
Using these, we proved the result: 
Theorem 2 Lei h(x) be analytic at xq, and 

f{x) = Wy lf(xo)l G ( '°°)- 



Then, 



h-\z)=X + f(x )J2^ n ~ 1 [f} [X 



,z - z 



I " 



ni 

71=1 



on a neighborhood of z = h(x ). 



In this paper, we study the asymptotic behavior of the nested derivatives 
£> n [/] (x) for large n, using a discrete version of the ray method [3], jl], [6]. 
As a consequence, we obtain asymptotic approximations for the higher-order 
derivatives of inverse functions. 



2 Nested derivatives 

Since ([T]) is difficult to study asymptotically, we obtained in [5] a linear 
relation between successive nested derivatives. 



2 



Proposition 3 Let 



9n[X) ~ [f(x)] n ■ { ) 



Then, g (x) = 1 and 

9n+i = g' n + in + 1) w(x)g n , n = 0, 1, . . . , (3) 

where 

= ff (4) 

As a result, we obtain the following corollary. 
Corollary 4 Let 

H(x) = h~ 1 (x), f(x) = ^-y, z = /i(:ro), |/(x )| G (0,oo) . 

T/ien, 

d n if 

"^T^o) = [/(^If^iK)- ra = l,2,.... (5) 

Later on, we will need to know the behavior of g n (x) for a fixed value of 
n and large x. 

Proposition 5 Suppose that 

ou(x) ~ ax p , x — > oo. (6) 

TTien, /or fixed n, we have 

r (-ir^i(-p-i)n ^ n+1 > p<-i 

<fe(*)~< (a-ir(^) n X- p=-l (7) 
[ n!a n a; pn , p > -1 

as x — >■ oo, where {p) n denotes the Pochhammer symbol defined by (p) = 1 
and 

ji-i 

(p)„=n(p + ^' «=i,2,.... 
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Proof. From ([6]), it follows that 

g n (x) ~ c n x Tn , x ->■ oo, 
for some sequences c n , r n . Since = u(x), we have 

ci = a, n = p. 

Using (J8| in ([3]), we get 

c n+ ix r,l+1 ~ c n r n x r ' l_1 + (n + 1) c n ax r,l+p , a; 

Comparing powers of x in (|9]), we obtain 

r n - 1, p < -1 
r n +p, p > -1 



Let's consider these cases one at the time. 



1. p < -1 
Solving 

we obtain 



r n+ i = r n - 1, n = p, 



r n = p + 1 — n. 



Using (10) in Mm, we must have 

c n+ i = {p + 1 - n) c n , ci = a, 



and therefore 

2. p > -1 
Solving 

we obtain 



c n = a(p + ir(-lT(-p-l) n . 

r n +i = r n +p, n= p, 
r n = pn. 



Using (11) in (|9j), we need 

c n+ i = (n + 1) ac n , a = a, 

and hence 



c„ = ma . 



3. p = -1 
We have r r , 



— n and 



which gives 



Cn+i = [-« + (n + 1) a] c n , ci = a, 



c n = (a - l) r 



a — 1 



Remark 6 Since 



71- 1 



" 3=0 

71-1 



fa-r 



a — 1 
[][a + (a-l)j], 



3=0 



we /lane 



(a - ir (^r). 



—7-1, as a 1. 



Hence, when p = — 1 and a = 1, we see that 

g n (x) ~ aT n . 



(12) 



In [7j we analyzed the family of polynomials generated by ^ with 



x. 



3 Asymptotic analysis of g n (x) 

We seek an approximate solution for (|3j) of the form 

G n (x) ~ /texp n) + n)] , n — >■ oo 
where k is a constant and 

G = o(F), n -> oo. 



(13) 



Since G (x) = 1, we require that 

F(x,0) = and G(x,0) = 0. 



Using (13) in <\3h, we have 



(f 9F ld2F G 9G } 
\ dn 2 dn 2 dn J 

dF dG\ 

— + — j exp (F + G) + (n + 1) u(x) exp (F + G) 



(14) 
(15) 



where we have used 



^ , v dF, , ld 2 F, , 
F(x,n + 1) = F(z,n) + — (x,n) + + 



From (15) we obtain, to leading order, the eikonal equation 
dF fdF\ 

_ + (n + 1Ml) _ exp (_) = o, 



(16) 



and 



ld 2 F dG 
2 dn 2 dn 



dG 

dx 



dF\ 
dn J 



exp 1 — I - 1 = 0, 



or, to leading order, the transport equation 



ld 2 F dG dG ( dF 
2 dn 2 dn dx ^ \ dn 



0. 



(17) 



3.1 The rays 



To solve (|16|), we use the method of characteristics, which we briefly review. 
:t order partial differential < 

$(x, n, F,p, q) = 0, with p 



Given the first order partial differential equation 

dF dF 
Q 



dx ' dn 

we search for a solution F(x,n) by solving the system of "characteristic 
equations" 

dx d$ dn d$ 
dt dp ' dt dq' 
dp d$ d$ dq d$ d$ 
dt 

dF d$ d$ 



dx ^ dF ' dt 



dn ^~dF' 



dt ^dp^^dq 
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with initial conditions 

3 [x(0, s), n(0, s), F(0, s),p{0, s), g(0, s)] = 0, 

and 

-f F(0, s) = p(0, s)-f x(0, s) + g(0, s)-fn(0, s), 
as as as 



(18) 



(19) 



where we now consider {x, n, F, p, q} to all be functions of the variables (t, s) . 

(20) 



For the eikonal equation (16), we have 



5 (x, n, F,p,q) = p — e q + (n + 1) u(x) 
and therefore the characteristic equations are 



dx dn 
~dt ~ ' ~dt 



dp , . I( v dq 

— [n + 1) oj (x), — = —oo[x) 



dt 



dt 



and 



dF 
~dt 



p — qe H 



(21) 
(22) 



Solving (21) subject to the initial conditions 

x(0,s) = s, n(0,s) = 0, q{0,s) = A{s), p{0,s) = B(s) 
with A(s),B(s) to be determined, we obtain 

x(t, s) =t + s, n(t, s) = exp [A(s)) f(s) [h(s) - h(t + s)) , 



p(t,s) = exp [A(s)\ 



lf(t + s) 
q(t,s) = In 



+ u(s) - (n + l)w(t + s) + 5(s), 
/(«) 



/(* + *) J 



+ A(s). 



From (18) we have B — e + w(s) = and therefore 



Thus, 



B = e A -co(s). 

x(t, s) =t + s, n(t, s) = exp [A(s)] /(s) [a(s) - h(t + s)) , 

/(«) 



p(t, s) = exp [A(s)] 
q(t,s) = In 



/(* + s) 

m 



f{t + *)\ 



-(n + l)w(t + s), 
+ A(s). 



(23) 
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Since ( 14 ) implies that 



F(0,s) = 0, 



we have from (19) and (23) 



[e A - u(s)] x 1 + A x = 0. 
Hence, A(s) = In [oo(s)] and therefore 

x = t + s, n = f'(s)[h(s) -h(t + s)} 



V 



/'(*) 



fit + s) 

3.2 The functions F and G 



n + l)u;(t + s), q = In 



/'(*) 



/(* + *) J 



Using (26) in (22) we have 



dF 



f'(s) 



dt f(t + s) 



(n + l)w(t + s) - In 



f(g) 



/'(*) 



/(* + *)' 



Solving (27) subject to (24), we obtain 



F(t,s) = In 



r m ' 


— n — n In 


|7(* + *)1 


[f(t + s)_ 







or, using (25), 



F = In 



/(f) 



n — n In 



/'(*) 



To solve the transport equation (17), we need to compute 
as functions of t, s. Use of the chain rule gives 



9" i ~n and 

on^ ' on 



' dx 


dx " 




-at 


at- 




"1 


0" 


$n 


ds 




dx 


ari 






dn 




ds 


ds 







1 


.at 


ds. 




.dx 


dn- 







and hence, 



-at 


m- 


1 


dn 


dx ' 


dx 


dn 








ds 
.dx 


ds 
dn- 


" J(t,s) 


§n 
. dt 


dx 
dt . 



where the Jacobian J(t, s) is defined by 

dx dn dx dn dn dn 



J(t,s) 



dt ds ds dt ds dt 
8 



Using (25), we find that 



J(t,s) 



/"(*) 



+ 0J(s). 



Using q 



|^ in (17), we have 



1 dq 
2d~n~ 



dG 
dn 



dG 



or 



d_ 
dn 



dx 

dG 

dx 



dG 
dn 



and using (21), we obtain 

d_ 
dn 



dG dx dG dn _ dG 
dx dt dn dt dt 



Since — e q 



we have 



d_ 
dn 



2 
1 

2J 



_ld_ 
2d~n~ 
d 2 n dx 

+ 



dn\ 
~dt) = 

d 2 n dx 



d 2 n dt d 2 n ds 
+ 



dt 2 dn 



dt 2 ds dtds dt 
1 d ( dn 
~2Jdt Ids ~~ 





1 / d 2 n 


)- 


2J V dt 2 


dn\ 


1 dJ 


~dt) 


~2J~di' 



dtds dn 
d 2 n \ 
+ dtds ) 



where we have used (30) and (31). Thus, 



dG 
~dt 



and therefore 



G(t,s) 



1 



\_dJ_ 
~2J~dl 

\n(J) + C(s) 



(32) 



for some function C(s). Since from (14) we have G(0, s) = 0, while (32) gives 
J(0, s) = co(s), we conclude that C(s) = | In [co(s)] and hence 



G{t,s) 



In 



[f'(s)f 



lf'(s)Y + nf(s)f»(s) 



(33) 
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Replacing (29) and (33) in (13), we obtain g n (x) ~ K$(x,n;s) as n — > oo, 
with 



$ (x, n: s) 



[f'(s)Y 

[fi(s)} 2 + nf(s)f"(s) 



and k is still to be determined. Eliminating t from ( 25 ) we get 

n - f{s) [h(s) - h{x)\ = 0, 

which defines the function s(x, n) implicitly. In cases where there exist mul- 
tiple solutions sx,S2,..., we must add all the contributions. 

We summarize the results of this section in the following theorem. 

Theorem 7 Let the functions g n {x) be defined by 

9n+i = g' n + (n + l)uj(x)g n , 

with go(x) = 1 and 



Then, we have 



= — In [f(x)] . 



g n {x) ~ K $ [x, n; Sj(x, n)] , n — > oo 



(34) 



where 



$ (x, n: s) 



m 

/(*) 



1 1 „\~\ n 



fix) 



[my 



[f'(s)Y + nf(s)f"(s) 



(35) 



k is an overall constant to be determined by matching and Sj (x, n) is a solu- 
tion of the equation 

n - f(s) [h(s) - h(x)} = 0. (36) 



4 Examples 

1. The natural logarithm. 
Let h(x) = ln(x + 1). Then, 



h'{x) 



x + 1 



(37) 
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and 

u(x) — (x + 1) 1 . 
In this case, ^ takes the form 

/ , n + 1 
9n+i = H TT#«' 9o = !■ 

X + 1 



(38) 



Using (37) in (35), we have 



$ (x, n- s) = (s + 1) e~ n (x + , 



while (36) gives 



/, - In ( | = 

x + 1 



or 



s = (x + l)e n -l. 



Thus, from (34), we obtain 



Since 



9n{ x ) ~ K ( x + 1) " ? Ti — )• OO. 



w(x) = (x + 1) 1 ~ x 1 , x — > oo, 



we know from (12) that gVi(x) ~ x n and therefore « = 1. We conclude 
that 

g n (x) ~ (x + l)~ n , n -)■ oo. 



But in fact, g n (x) = (x + 1) ™ is the exact solution of (38)! 

2. The arctangent. 

Let h(x) = arctan(x). Then, 

f(x) = x 2 + 1 

and 

u)(x) 



(39) 



2x 



In this case, ^ takes the form 



x 2 + 1 



2x 



9n+l = 9' n + (n + 1) 2 . . 9n, 90 = 1- 
X z + 1 
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Using (39) in (35), we have 

s 2 + l 



$ (x, n: s) 



x 2 + l 



2s 



x 2 + l 



2s 2 



2s 2 + n(s 2 + 1) 



(40) 



while (36) gives 

n — 2s [arctan(s) — arctan(x)] = 0. (41) 
For every point (x, n) there exist two solutions s_ < and s + > of 



(41) (see Figure 1). Hence, we get from p4| ) 

g n (x) ~ k [$ (x, n; S-) + $ (x, n; s+)] , n — > oo. 



If we fix n = 5 and let x — > oo, we get from (41 ) 



X 


s 


10 


-1.0877, 35.114 


20 


-1.0686, 70.057 


50 


-1.0575, 175.02 


100 


-1.0538, 350.01 



It follows that one solution approaches a fixed negative value s_ and 
the other s + grows algebraically. After some calculations, we find that 



s + = (l + D x + O (x x ) , x -)■ oo. 



Using (42) in (40) we obtain, to leading order, 



$ (x, n; s+) ~ k2~5 (n + 2) n+ * e~ n x' n , x ->■ oo. 



(42) 



(43) 



But since 



we have from ([7]) 



cj(x) = — r-^ — ~ 2x 1 , x — > oo 
x" 1 + 1 



<?n (#) ~ (^ + ™, x — > oo. 



(44) 



Matching (43) with (44), we get 



K = 22 (n + 2y n ~*e n (n + l)\ 



12 
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Figure 1: A sketch of s-(x,n) and s + (x,n) for various values of n. 
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Note that 

K = 4 v / 7re -2 + O (n^ 1 ) , n -)■ oo. 

We conclude that 

SW^e) ~ 4A/7re~ 2 [$ (x, n; s_) + $ (x, n; s+)] , n — >■ oo. 



If x = 0, (41) becomes 



n — 2s arctan(s) = 0, 



(45) 



and solving for s, we get 



n 


s 


5 


±2.1887 


10 


±3.8056 


20 


±6.9985 


50 


±16.551 


100 


±32.467 



It follows that in this case s_ = — s+, and therefore 
n (O) ~ 4v^e- 2 [$ (0, n; s+) + $ (0, n; -s+)] 

= 4v^[l + (-1)™] (4 + 1) e-(" +2 ) (2s+) n 



2s 2 , 



2s 2 ±n(4±l) 



as n — > oo. Since g n (0) = for odd n, we focus our attention on even 
values of n. From ( [45] ) we obtain 

n + 2 

Si ~ , n — y oo, 



7T 



and hence, 



^2n(0) 



4 2(n+l) 
7T 2 "+f ' 



-2(n+l) 



(n + 1) 



2n+1 4(n + l) 2 + 7r 2 



n — )■ oo. 



4(n + l) 3 + 7T 2 ra 



But since in this case if (x) 
£ 2 "[x 2 + 1] (0) 



tan(x), we know that [2] 
2 



n + 
14 



_ 4 - ( 4 n+i _ !) | Bz(n 



+i)| > 



where B n are the Bernoulli numbers. It follows that 

4« +1 f n + l\ 2{n+l) 4(n+l) 2 + 7r 2 



|-B 2 ( n+ i)| 



^ + l(4^ -1)V e J v / 4(n + 1) 3 + 7r2n 

as n — > oo, or 

.„ n 4 n fn\ 2n An 2 + -K 2 

Bo n \ ~ 2— — = - : n ->• oo. 



5(4 n -l)Ve7 ^n 3 + TT 2 (n - 1)' 



To leading order in n, we obtain the well known asymptotic approxi- 
mation 

|S 2n | ~ 4v / nvr ( — ) , n — )> oo. 
Ve7r/ 
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